typeset using JPSJ.sty <ver.l.Ob> 



Effects of Orbital Degeneracy and Electron Correlation 
on Charge Dynamics in Perovskite Manganese Oxides 



On 
ON 
On 



ON 



c3 



Hiroki Nakano, Yukitoshi Motome 1 and Masatoshi Imada 

Institute for Solid State Physics, University of Tokyo, Roppongi 7-22-1, Minato-ku, Tokyo 106-8666 
1 Department of Physics, Tokyo Institute of Technology, Oh-okayama 2-12-1, Meguro-ku, Tokyo 152-8551 

(Received ) 

Taking the orbital degeneracy of e g conduction bands and the Coulomb interaction into account 
in a double-exchange model, we investigate charge dynamics of perovskite Mn oxides by the 
Lanczos diagonalization method. In the metallic phase near the Mott insulator, it is found that 
the optical conductivity for a spin-polarized two-dimensional system exhibits a weight transfer 
to a broad and incoherent structure within the lower-Hubbard band together with a suppressed 
Drude weight. It reproduces qualitative feature of the experimental results. As an orbital effect, 
we find that an anomalous charge correlation at quarter filling suppresses the coherent charge 
dynamics and signals precursor to the charge ordering. 
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Although properties of perovskite manganese oxides 
i?i„ x A x Mn0 3 (i?=rare earth, A==Ca, Sr, Ba or Pb) 
showing colossal magnetoresistanceEP are to some extent 
reproduced by a double-exchange (DE) model composed 
of non-interacting (NI) electrons in a non-degenerate e g 
conduction band ferromagnfitically coupled with t 2 g lo- 
calized spins of 5=3/2 JrElertm there remain open prob- 
lems which cannot be understood within the framework 
of this model. One of them is that Lai_ x Sr x Mn03 in 
the ferromagnetic and metallic phase near the insulator 
(:c=0.175) shows a large and broad structure of incoher- 
ence in the optical response together with a relatively 
small Drude weight at low temper atures0EF though the 
spin degrees of freedom cannot contribute to the inco- 
herence under the perfect polarization of spins. Some 
other degrees of freedom not contained in the DE model 
must be involved. Possible origins of the incoherence are 
the dynamical Jahn- Teller (JT) distortion,™ ^orbital 
degeneracyEj'OE3 ) and electronic correlations In 
ref. |^ employing the dynamical mean-field approxima- 
tion, the electronic correlations are not taken into ac- 
count though they are not negligible and play an essen- 
tial role as shown later. In this work, we focus on the 
latter two and make it clear how large incoherence is 
induced in a minimal model of these two origins. 

On this subject, related studies have already been done 
under some approximations. Shiba, Shiina and Taka- 
hashi claimed that the large incoherence in the optical 
absorption originates in the degeneracy of the e g or- 
bitals,El although the Mott insulator in the undoped 
system is not reproduced becausG-tha,s±xong correlation 
effect is not seriously considered J19> t3c3 Ishihara, Ya- 
manaka and Nagaosa derived an effective t-J type Hamil- 
tonian considering thej-degeneracy of the e g orbitals and 
the strong interaction.!!^ They obtained incoherent part 
of the optical conductivity in a mean-field treatment. In 



such a mean-field treatment, unfortunately, the Drude 
weight is absent with a diverging specific heat coefficient 
7 in contrast to the experimental resultsH Horsch et 
al. applied a finite-temperature Lanczos method to the 
'orbital' t-J modelxfP In their method, however, the 
accessible temperature range is much higher than the 
experimental conditions JHP Consequently their Drude 
weight is largely suppressed simply by thermal fluctua- 
tions, which makes hard to compare with the experimen- 
tal results. In order to discuss this metal-insulator tran- 
sition, fluctuation effects are important; a reliable way 
to calculate beyond biased approximations is required. 
In this work, we examine the Hamiltonian given by 



"H=V V t vh 



cJ„Gj 



(1) 



as a minimal model mentioned above. This model is 
derived from the generalized DE model with e g -orbital 
degeneracy under the assumptions of strong Hund's-rule 
coupling and perfect spin polarizationO Here tJJ de- 
notes the hopping integral and U is the effective interor- 
bital Coulomb interaction obtained after subtracting the 
Hund's-rule coupling energy between e g electrons. The 
orbitals d x 2_ y 2 and d^ z i_ r 2 correspond to v—\ and 
v = 1, respectively. Note that this model is essentially 
different from the usual Hubbard model. Actually, the 
difference appears even in the NI case as shown later. Wc 
consider here only the nearest-neighbor hopping given 
by t}} = -3/4t 0) tj] = -l/4t , = -(+)V3/4t 

along x(y)-direction in two dimensions. Here, we inves- 
tigate the two-dimensional (2D) model given by eq. (|l|) 
based on the following reasons. One is that the three- 
dimensional (3D) system of R\- x A x yhiO-i keeps a per- 
fectly spin-polarized plane not only in the ferromagnetic 
doped phase but also in the ^4-type antiferromagnetic 
insulating phase at x—0 so that the spin polarization is 
retained in a 2D plane over all regions of interest. The 
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second is that a similar incoherent feature in the optical 
conductivity was also reported in the 2D system, for in- 
stance, in La1.2Sr1.gMn2O7.li3' It is therefore likely that 
the 2D system can capture basic physics of the incoher- 
ent charge dynamics. The model ([!]) contains an implicit 
variable parameter of a doping concentration 5 = 1 — & 
or an electron density n = where N c and N s de- 
note the numbers of electrons and sites, respectively. 
Recently, this model (Q) has been studied by—means of 
quantum Monte Carlo (QMC) method.BBE]) There, 
the relation was discussed between the incoherence of the 
charge dynamics and divergences of both a staggered or- 
bital correlation length and the compressibility as 6 — > 0. 
To get more insight, it is desired to calculate quantities 
directly connected to the charge dynamics. In this pa- 
per, then, we discuss effects of the orbital degeneracy and 
the electronic correlation based on direct calculations of 
dynamical quantities about charge transport. 

We here employ, as another unbiased method, ex- 
act diagonalization of finite-size (FS) clusters, using the 
Lancz6s_algorithm and the continued-fraction-expansion 
methodEr to obtain dynamical quantities at zero tem- 
perature. The sizes of the cluster calculated in this work 
are \/T0 x yTO and 4 x 4-Jn the y/lO x vTO-site system, 
to reduce the FS effect,Ef we introduce the Aharonov- 
Bohm flux $ in the hopping term as the twisted bound- 
ary condition and find the optimized boundary so as to 
minimize the ground-state energy. In the 4x4 system, 
we choose the boundary condition among the periodic 
($ = (0,0)), the anti-periodic (<& = (tt,tt)) and the 
mixed (<& = (0,tt) or (ir, 0)) ones. 

We here calculate the optical conductivity defined as 
a(u)=[a x (uj)+ay(u!)}/2 where cr a (u>) = 2ite 2 D a 5{uj) + 

TrE m(/ o) K Slr^ w ~E m + E Q ). Here, D a is the 
Drude weight; j a is a current operator along a-direction 
(a=x,y) given by j a = -i E ^,^5 J c L c i+<5 Q ,^' ~ 
cj , g v ,Civ); and \m) denotes an eigenstate of the sysytem 
with the energy eigenvalue of E m . Note that m = 
represents the ground state. The averaging operation 
is performed to reduce the FS effect and to cancel the 
anisotropy due to the anisotropic boundary condition. 
The sum rule / °° a(uj)dui — ne 2 K is satisfied, where 
— AK is the kinetic energy per site. We will call K 
the total weight hereafter. When U is large enough, 
we also calculate the effective carrier density defined as 
N e g — ^2 Jq c (j(uj)duj where cj c is a frequency just be- 
low the bottom edge of the upper-Hubbard (UH) band. 
In this work, we make a further procedure to reduce FS 
effects. The procedure is a multiplication of a correc- 
tion factor determined from the FS effects in the NI 
cases. The factor is defined by r = K^~^ / K^~°\ 
where K^~ ^ and K^~ ^ are the total weights in the 
thermodynamic limit and in the 7V s -site case when U=0, 
respectively. Hereafter, all the quantities obtained af- 
ter this procedure are labeled with suffix c, for example, 
D c . This procedure successfully reduces the FS effects-in 
the usual Hubbard model in one and two dimensions z3 
Here, we take a large value of U/to—16, at which a large 
Hubbard gap opens at half filling and the FS effect is 
expected to be less serious. 



We first show our results of the doping dependence of 
the total weights, the effective carrier densities and the 
Drude weights in Fig. 1. Even in the NI case, as was re- 
ported in the 3D model in ref. |l(], the present model ([!]) 
in 2D exhibits a different doping dependence between the 
Drude weight and the total weight due to the hopping 
between the d x 2_ y 2 and d 3z 2_ r 2 bands. (See the gray 
curves in Fig. 1.) However, the Drude weight remains 
finite even at half filling, which indicates that the system 
is metallic in this NI case. When U/to = 16, on the other 
hand, the Drude weight at half filling has a very small 
value ~ 0.0005 in the case of vTO xVTO-site system. This 
small value indicates that the system is Mott insulating 
and that even a cluster of VTO x \/T0 sites can reproduce 
well the bulk quantities. In the dilute-electron-density 
region (8 ~ I), the total weights, the effective carrier 
density and the Drude weights for U/to — 16 behave 
close to the corresponding quantities for the NI case, 
which indicates that effects of Coulomb interaction are 
small in this region. At 6 — 0.5, the Drude weight shows 
a significant dip, which reminds us of a charge ordering 
with a Ci?-type magnetic structure observed in such-ma- 
terials as Nd .5Sro. 5 Mn0 3 E|) and Pro.5Cao.5MnO3.BE3) 
Although the value of the Drude weight is still finite, this 
dip is neither seen in the NI case of the present model 
nor in the usual Hubbard model with finite Ux3 We will 
discuss this issue later. With decreasing S further, the 
Drude weight and the effective carrier density vanish for 
U/to=16 while the two quantities increase for U=0. It 
should be emphasized here that the overall behaviors of 
the Drude weight and the effective carrier density except 
for the dip at 6=0.5 in the Drude weight are qualita- 
tively similar to the ones of the usual Hubbard model in 
2D.E3 We, however, also note that the critical region of 
suppressed Drude weight near 6=0 seems to be narrower 
than that in the 2D Hubbard model due to the imbal- 
anced populations of electrons in d x 2_ y 2 and d 3z 2_ r 2. 

In Fig. 2, we present results for the incoherent part of 
<j{u>). One can see in Fig. 2 (a) that, at half filling, a 
large gap exists and that only the weight transfer across 
the gap to the UH band appears. As shown in Figs. 2 
(b) and (c), the more holes from half filling are doped, 
the more weights are transferred from the region above 
the gap to the region within the lower-Hubbard (LH) 
band. Especially, note that weights for [//io=16 spread 
in a wider region than those in the NI case shown in the 
inset of Fig. 2(b). The QMC study shows a critical en- 
hancement of the orbital correlation length induced by 
U whictuwould cause the incoherence of the charge dy- 
namics.EZP The orbital correlation in the present calcula- 
tion is in agreement with that in the QMC. The present 
calculations indeed show a close relation between the or- 
bital fluctuation and the mid-gap incoherence. The qual- 
itative feature of the incoherent charge response exper- 
imentally observed in <j(ui) is reproduced in the present 
calculation, which supports the importance of combined 
effects from orbital degeneracy and strong correlations. 
On a quantitative level, however, D/K ^0.51 at 6 ~0.125 
in the present result is still larger than the experimental 
indications, ~0.2. 

To understand the dip at (5 = 0.5 in Fig. 1, let us dis- 
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cuss the orbital_and charge structures in the ground 
state at S = 0.5.E3 The basic point is that an orbital-po- 
larized state with a spatially anisotropic overlap seems 
to enhance the charge ordering as we see below. To 
gain the kinetic energy in such a state at finite U and 
S = 0.5, an occupied site favors a vacant site as the near- 
est neighbors. In the model (Q), the orbital polariza- 
tion (OP) defined by ((n x 2_ y 2) - (n 3z 2_ r2 ))/((n x 2_ y 2) + 
(n 3z 2_ r 2)) grows because d x 2_ y 2 orbital has larger hop- 
ping matrix than d% z 2_ r 2 orbital— .The strong interaction 
also enhances the polarizationO Actually the OP is 
- 0.756 at (5 = 0.5 for J7/t =16E§ in the 4x4-site sys- 
tem. Note here that, on the other hand, the usual 
Hubbard model does not show such a spin polarization 
in the ground state at (5 = 0.5. When one uses \0) = 
cos 6\d x 2_ y 2) + sin 9\d 3z 2_ r 2) as a notation for a single- 
site state, the above polarization leads to \0\/tt ~ 0.114, 
where each occupied orbital is strongly anisotropic and 
has a large overlap only in the x- or y-direction. We have 
next calculated the charge density correlation defined 
as A(fe) =N- 1 Y, i , j [(n i n j ) - (ni)(rij)] exp[ifc • ( Xi - Xj)] 
in the 4x4 sites, where n% =nn +n;2- The results at 
S = 0.5 show a peak at fc = (7r,7r) which makes stag- 
gered charge correlation (CC). The doping dependence 
of the above CC at (tt, it) is shown in Fig. 3. An anoma- 
lous peak appears at (5 = 0.5 only for the model ([!]) while 
the usual Hubbard model exhibits a monotonic 5 depen- 
dence. A similar behavior at (5 = 0.5 of the CC and the 
charge dynamics was reported in ref. where a spinless 
and orbitless system with the usual Hubbard-type hop- 
pings and the nearest-neighbor interaction V was stud- 
ied. Although the size treated here is too small to judge 
whether the present system is metallic or insulating, an 
anomalous behavior of A(tt,tt) at (5 = 0.5 occurs even in 
the present case without V. The enhanced CC at 6 = 0.5 
which makes the suppressed coherence in the charge dy- 
namics can be induced only by the anisotropic hopping 
and the on-site interaction while the static charge or- 
der will be more stabilized with the help of the inter- 
site Coulomb replusion V. The JT distortion, stud- 
ied theoretically by means of LDA+t/o- 1 and Hartree- 
FockM methods, will also enhance the above tendency 
of the CC and the incoherence of the charge dynam- 
ics. Although three dimensionality would reduce the OP, 
which influences the stability of this staggered CC, our 
results in 2D system are important because the charge- 
ordering phenomena essenti ally . a npear with strong 2D 
anisotropy in experiments.E3ESE3'E2|) To obtain infor- 
mation of orbital patterns at (5=0.5, in addition, we cal- 
culate orbital correlations defined by C = ^^(WiWj), 
where Wi denotes a psuedo-spin operator of 3x 2 — r 2 
or 3y 2 - r 2 {Wi = -%T? + (-)^Tf for 3x 2 (y 2 ) - r 2 , 
and Tf = ± J2 VV < Ku' c iv c iv' with the Pauli matrix 6% v ,). 
Prime at the sum means that i and j run over one sub- 
lattice. The possible patterns on the lattice are illus- 
trated in Fig. 3. The results reveal that C for (a) is 
the largest£3) The_pattern (a) agrees with the exper- 
imental indication. lj' Thus, the basic orbital structure 
in charge-ordering phenomena is determined irrespective 
of magnetic structure due to difference of energy scales. 



This situation is captured well within the model (|Tj) in 
which the spin degrees of freedom are frozen out. 

Finally, it is worth mentioning effects of the three di- 
mensionality neglected above. In 2D, there happens the 
imbalance of populations in d x 2_ y 2 and d 3z 2_ r 2 as de- 
scribed above. In 3D case, a hopping along z-axis would 
reduce the imbalance. We note that, in the NI case, such 
a reduction of OP in 3D from that in 2D decreases the 
coherence as inferred from the comparison between 3D 
case in ref. [l^ and 2D one in Fig. 1. This is simply be- 
cause the charge incoherence can not be induced if the 
orbital is completely polarized. The orbital depolariza- 
tion may also be simulated by introducing a chemical 
potential difference between the two orbitals to compen- 
sate the existing OP. We have performed calculations for 
the VlO xy/lO -site 2D Hamiltonian with such a chemi- 
cal potential difference. The results show that the Drude 
weight becomes more suppressed. For example, if the 
chemical potential difference is tuned to keep the vanish- 
ing OP, ~ 0.45 at 5=0.2 as compared to ~ 0.67 
in the case of Fig. 1. The broad shape of incoherence 
within the LH band near the Mott transition is found to 
be qualitatively unchanged after the chemical potential 
is introduced. These suggest the possibility that the de- 
polarization due to the three dimensionality would also 
make the charge transport more incoherent. 

In summary, we have investigated the optical conduc- 
tivity (j(ui) of the double-exchange model with both the 
e s -orbital degeneracy and the Coulomb interaction. As- 
suming the spin polarization and the strong Hund's-rule 
coupling in the 2D system, we have examined the effects 
of the orbital degeneracy and the on-site interaction. We 
have obtained in cr(u>) a suppressed Drude weight and a 
mid-gap incoherence with a broad structure in the metal- 
lic state near the Mott transition. The results reproduce 
to a considerable extent the incoherent charge dynamics 
observed experimentally in Mn oxides. On a quantitative 
level, however, the calculatecLDrude weight does not ap- 
pear to be sufficiently small. lIlP Effects of fluctuations 
of JT distortion would be an additional source of the 
incoherence, which is an interesting subject of further 
studies for more quantitative comparison. We have also 
found an anomalous feature in charge transport and the 
charge correlation at quarter filling, which reproduces 
basic structure of the orbital and charge ordering ob- 
served in experiments. 
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Fig. 3. Doping dependences of the charge correlation. An 
anomaly at 5=0.5 for the model (jy) is obtained while only the 
monotonic behavior is seen in the usual Hubbard model. Insets 
(a), (b) and (c) display possible staggered patterns of orbitals 
at 5=0.5 in 4x4 cluster. Simple circles and anisotropic symbols 
denote Mn 4 + and Mn 3+ with 3x 2 — r 2 /3y 2 — r 2 , respectively. 



Fig. 1. Doping dependences of (a) the total weight K, (b) the ef- 
fective carrier density N c ff and the Drude weight D at [7/trj=16. 
For comparison, gray lines show results in the non-interacting 
case, where the total weight, the effective carrier density and the 
Drude weight in the thermodynamic limit jV s — > oo are shown by 
the dotted line in (a), dotted one in (b) and solid one in (b), 
respectively. 



Fig. 2. Incoherent part of the optical conductivity at U/to = 16 
for (a) the half-filled case in the 10-site system under the anti- 
periodic BC, (b) 2 holes in the 4x4 sites under the mixed BC, 
(c) 4 holes in the 4x4 sites under the anti-periodic BC. Inset in 
(b) shows the non-interacting case at the same filling and the 
same BC. Delta functions are broadened with width of 0.05<o. 
We choose the BC to realize the lowest-energy ground state. 
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